Liquid-like spatial distribution of magnetic polarons revealed by neutron scattering in 

La^Ca^MnOs 
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Elastic neutron scattering experiments performed in semi-conducting Lai-zCa^MnOs single crys- 
tals (x=0.05, 0.08), reveal new features in the problem of electronic phase separation and metal 
insulator transition. Below Tat, the observation of a broad magnetic modulation in the q-dependent 
scattering intensity, centered at nearly identical q m whatever the q direction, can be explained by 
a liquid-like spatial distribution of similar magnetic droplets. A semi-quantitative description of 
their magnetic state, diameter, and average distance , can be done using a two-phase model. Such 
a picture can explain the anomalous characteristics of the spin wave branches and may result from 
unmixing forces between charge carriers predicted from the s-d model. 
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Doped Mn perovskites are intensively studied for their 
remarkable giant magnetoresistance underlying the elec- 
tronic behaviour that has potential technological applica- 
tions. The basic physical ideas have been given by Zener 
HI, with the model of the double exchange, which con- 
nects the electron hopping with the ferromagnetic (F) 
alignments of Mn spins. The electron-phonon coupling 
through the Jahn- Teller effect was studied by Millis et 
al. On the other hand, the role of the electronic lo- 
calization and the formation of magnetic polarons in this 
problem has been emphasized by Varma et al. ]3|. Sev- 
eral pictures of magnetic polarons were proposed, associ- 
ated to one carrier (small polarons) |4|,|| or, in the case 
of larger carrier densities, to several carriers (large po- 
larons) ||. A cooperative state of magnetic droplets by 
self-trapping carriers was predicted by Nagaev using the 
s-d or s-f model , since this physical situation was first 
suggested for some rare earth compounds . Models of 
electronic phase diagram have recently been reexamined 
HP]. Magnetic inhomogeneities have been suggested to 
be at the origin of the behaviour of the susceptibility and 
of the dynamical spin fluctuations observed above T c in 
doped Mn perovskites jR)|,|ll| . Very recent NMR experi- 
ments have been interpreted in terms of electronic phase 
separation However, a characterization of magnetic 
inhomogeneities has never been made until now. We have 
started a spin dynamics study in Lai-xCa^MnOa for low 
doping Jl4|,[l5| where the system is insulating and under- 
goes a weakly canted antiferromagnetic (AF) transition 
at Tn [ p"3[ . An additional spin wave branch occurs by 
doping. The corresponding dynamical susceptibility re- 
veals a small ferromagnetic correlation length, whereas 
the dispersion curve is isotropic, indicating a new ferro- 
magnetic coupling between some Mn spins. This isotropy 
contrasts with the strongly anisotropic character of the 
other spin wave branch, characteristic of super-exchange. 
This new spin dynamics was attributed to "magnetic po- 
larons" , whose origin could not be elucidated. 



We have performed elastic neutron scattering measure- 
ments as a function of temperature in Lai-^Ca^MnOa 
(x=0.05, 0.08), close to the direct beam (t = 0) and to 
t=(110) Bragg peak. In both experimental broad 
modulation is observed, centered at nearly the same q m 
whatever the q direction, with intensity growing below 
Tat. This scattering does not exist for x=0. Such a 
pattern is typical of an assembly of magnetic clusters 
or droplets, with a mean magnetization different from 
that of the matrix and a well-defined shortest distance 
between them. A semi-quantitative description of the 
droplets and of their spatial distribution is proposed us- 
ing a liquid-like model. The droplets are isotropic, their 
radius is w 9 A, their density is low compared to that of 
the hole concentration and their minimal distance of ap- 
proach indicates a repulsive interaction. The existence of 
these magnetic clusters is closely linked with the unusual 
spin dynamics (T^hTtJ]. The overall static and dynamic 
observations suggest a picture of magnetic droplets or 
large polarons with a magnetic coupling distinct from ex- 
change, coupled through the surrounding medium char- 
acterized by the super-exchange coupling. They are char- 
acteristic of an electronic phase segregation, which results 
from the balance between unmixing forces, predicted by 
the s-d model, and intermixing Coulomb forces. 

Experiments have been carried out at the reactor 
Orphee (Laboratoire Leon Brillouin) using a triple axis 
spectrometer set on a cold neutron source. The sam- 
ples with x=0.05 and 0.08 show similar results, so that 
only observations for the x=0.08 sample, with a volume 
of 0.3 cm 3 are reported. The crystalline structure is 
orthorhombic (O') with Pbnm symmetry. The crystals 
are twinned basically having three domains of volume 
Vj (£=1,2,3) with their c axes along the directions ao, 
bo, Co of the cubic perovskite substructure (called F do- 
mains in Jl4[). The mirror or M domains, leading to two 
peaks in the rocking curve, may have some influence close 
to t ^ 0. The q directions [110], [001], [112] and [100] 
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FIG. 1. a) Scattering Intensities versus q, calibrated in 
barns, observed above Tat (I(q)=I(q)Hr), at 75K and 15K for 
[110]i q direction, b) Magnetic intensities I m (q)=I(q)-IirT(q) 
versus q at 15K for [110]i, and [112]i q directions. The 
dashed lines are calculated functions according to the model 
described in the text. 
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FIG. 2. Diffuse scattering intensities versus Q (or q=Q-r) 
measured around r=(110)i. a) along [110]i, at five tempera- 
tures in the 300K-17K range (Tjv=122K). When not shown, 
the error bar is within the size of the symbol. In the inset, 
the continuous lines are components of a fit with two (+q m , 
-q m ) modulations, symmetric with respect to (110)i. b) At 
15K along three q directions. The horizontal line locates the 
sample background observed at 300K (cf fig 2-a). 



refers to a , c , to the diagonals of the faces (a +c ) 
and (ao+ bo) respectively. The subscript i will be used 
below to distinguish the q direction within the different 
F domains. The magnetic structure has been determined 
from Bragg peak intensities. It consists of ferromagnetic 
(a-b) planes, with AF stacking along c. However, the in- 
crease of the (110) and (112) nuclear Bragg peaks below 
T^ (122K), indicates a long range ferromagnetic order- 
ing for a small spin component along c. It corresponds to 
a canted state. At 14K, the average canting angle, 9 av , 
is «10 degrees from b within the (b-c) plane. 

Elastic experiments have been performed around r = 
with k,=1.25A , covering 6 temperatures within the 
15K< T <300K and 0.05< q <0.6l-l ranges, along 
several q directions. Intensities have been put on an ab- 
solute scale, using a vanadium sample and a standard 
procedure for correction. The temperature dependence 
of the spectra are reported in Fig 1-a for q/ /[110]i. In the 
150K < T < 300K range, I(q)= I(q) ht is found temper- 
ature independent. In the smallest q range, I(q)ijr is at- 
tributed to dislocations or large structural defects present 
in the sample. At larger q, the residual intensity is nearly 
twice as larger as the nuclear and chemical incoherent 
scattering (cf the arrow in Fig 1-a), indicating a contri- 
bution from disordered spins, static at the experimental 
time scale (0.5 10~ n s). Below Tjy, I(q) increases in the 
whole q range, keeping values very close to I(q)/rr around 
qwO.OlA -1 , so that a modulation appears centered at 
q m «0.2A _1 at 15K. The temperature dependence of q m 
cannot be detected in the 75-15K range within our ac- 
curacy. Similar observations are obtained along [112]i 
(mixed with [100] 3), except for a small shift of q m towards 
a smaller value and a small overall decrease of intensity. 
I(q) along [001]i (mixed with [110)2,3), is identical to I(q) 
along [1 10] 1 . The same observations are obtained using 
an XY detector device (SANS), at 15K where the dy- 
namical spin fluctuations can be neglected. In addition, 
in the smallest q range (q<0.08A _1 ) these SANS inten- 
sities show an enhancement along [112]i, which is typical 
of (112) interface planes between twin domains. 

Fig 1-b reports the magnetic contribution I m (q) 
obtained by substracting I(q)jjy, for q//[110]i and 
q//[112]i at T=15K. The smallest q range correspond- 
ing to effects related to the dislocations, is not shown. 
Im(q) consists of a broad peak, located at q m « 0.21A -1 
and a residual and nearly flat scattering beyond 0.4A -1 . 
Along [1 12] 1 , the overall intensity is slightly smaller, and 
a slight shift of q m is observed. 

Elastic experiments (ki — 1.55A -1 with Be filter), 
have been also performed in the ([110], [001]) plane close 
to r=(110)i, at 8 temperatures in the 11K < T < 300K 
range, and several q directions. The [110]i q direction 
is reported in Fig 2-a. In this figure, the two sets of 
vertical lines are results of the fit of the (110)i (bold) 
and (002)2 Bragg peaks, easily resolved due to the large 
orthorhombicity. At 15K, a well defined maximum in 
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I(q) is observed at C=0.87 (i. d. q m =0.22 _1 ), with a 
fiat q scattering beyond. The q-symmetric intensity of 
this modulation with respect to (110)i is observed as a 
shoulder on the huge (002) 2 Bragg peak side of the twin 
domain, as shown in the inset of Fig 2-a. The apparent 
asymmetry of the spectrum, is due to the superposition 
of (002)2 and (H0)i Bragg peaks with no modulation 
from the domain (002)2- We conclude that this modu- 
lation is magnetic, with a magnetization along (001). In 
Fig 2-b, I(q) is reported for q along two other directions 
of the scattering plane at 15K, taking the same q origin 
at (110)i. A small overall decrease of the intensity and 
a shift of q m to smaller values is observed. Along [1 12] 1 , 
the q m value appears slightly smaller than that observed 
close to r = ( Fig 1). This is well explained by the 
shift of the q origin attached to the different domains, 
unlike the previous case, close to r = 0, where the q ori- 
gin is the same for all domains. The intensity is smaller 
by a nearly q independent value. Actually, the intensity 
of the flat q tail keeps the low value determined at 300K 
(cf the horizontal lines in Fig 2 (a) and (b)) whereas 
along the other directions it varies with temperature. As 
the temperature increases, a very slight decrease of q m 
can be detected whereas the intensity of the modulation 
decreases (Fig 2-a). At T^, a huge critical scattering 
intensity is observed close to the (110)i Bragg peak. In 
the Tn< T < 250K range, small modulations appear, 
but along [1 10] 1 (ao axis) only. As shown at 250K in 
Fig 2-a, the q maxima correspond to the harmonics of 
q m . Their intensities show irreversibilities with temper- 
ature, in contrast with the broad modulation below Tjv 
which is perfectly reproducible. Since undetected close 
to r=0, they are likely related to modulated strain effects 
(with a Q=r+q dependent intensity) coexisting or com- 
peting with the broad magnetic modulation. At 300K, 
I(q) is monotonous, but its weak q dependence indicates 
the persistence of ferromagnetic correlations, static at 
the experimental time scale. These observations confirm 
those obtained close to the direct beam. The flat q scat- 
tering, with minima of intensity along [112], appears as a 
"sample background" corresponding to disordered spins. 
Although not explained, the minima of intensity along 
[112] could be a complementary effect of the maxima ob- 
served by SANS for q < 0.lA _1 , and attributed to (112) 
interfaces between domains. 

We focus our analysis on the magnetic contribution 
Im(q) obtained close to r = at 15K (Fig 1-b). A scat- 
tering pattern showing a broad modulation with nearly 
the same q m in all q directions, indicates a spatial or- 
ganization of similar entities or " droplets" , as that ob- 
served e. g. in a chemical unmixing process (lSy . To 
apply this model in a magnetic system, needs to define 
a mean magnetization density into the two phases, with 
a "contrast" between them, as an analog of the chem- 
ical contrast. The existence of a ferromagnetic compo- 
nent along c, allows to define a magnetization function 



m(r), m||c throughout the lattice. Therefore, the spatial 
Fourier transform of the spin correlations < SiSj > is re- 
placed by that of < m(r)m(r') >, r being a continuous 
variable in direct space. Neglecting the small anisotropy, 
a semi-quantitative analysis can be made using: 
da(q)/dn = rlA.N v .Vl\Am\ 2 .\F( q R) 2 \.J(q) (1) 
where Am is the difference between the mean magneti- 
zations m||c inside and outside the droplet. F(qR) is the 
form factor of the droplet (a spherical shape of radius R 
is assumed) of density Ny and volume V<; and J(q) is the 
interference function. The A factor (A=T,iViSin 2 9i/'EiVi 
where 9i is the angle between the mean magnetization m 
and q for the domain V^) is required to determine Am. 
It expresses that for a q direction belonging to domain 1, 
the intensity also comprises contributions from domains 
2 and 3. They have different q and m directions and are 
weighted by their geometrical factor sin 2 8i. This latter 
factor shows that [110] is a "pure" direction, [001] can- 
not be observed and [112] is mixed with [100], weighted 
by 0.5(Vi+V 2 ) and V3 respectively. The equality of the 
three twinned volumes could be checked at 11K from the 
scattering pattern of the SANS experiment, thanks to 
the XY detectors device. We have used the isotropic 
J(q) function derived for liquids by Ashcroft and Lekner 
Jl9| . The spatial distribution is characterized by two pa- 
rameters: the particle density Ny from which we deter- 
mine the mean interparticle distance d m and the minimal 
distance of approach between the particle centers d m i n - 
The expression of J(g) is given in the footnote pp| . This 
model can fit the experimental data as shown in Fig 1- 
b (dashed line). It yields a diameter 2R « 18 A (i.e. 
«4-5 lattice parameter ao of the small perovskite cube) , 
a droplet density Ny = 2.1 x 10~ 5 A~ 3 leading to d m 
« 9 ao, and a minimal distance of approach d m i n = 26 A 
(6 ~ 7ao). This value, compared to the average dis- 
tance d m , accounts for the rather well defined organi- 
zation and indicates a repulsive interaction between the 
droplets. From the absolute value of the intensity, we 
determine a magnetic "contrast" Am 0.6 ± 0.20/is. It 
corresponds to a difference of «8 degrees, between the 
canting angles characteristic of the two phases, which is 
surprisingly small. 

From this semi-quantitative determination, we con- 
clude that the density of the "ferromagnetic" particles 
is very low compared to that of the holes (ratio 1 /60 for 
x=0.08) leading to a picture of hole-rich droplets within a 
hole-poor medium B||E|- At x=0.08, the magnetic state 
of the droplet is far from a true ferromagnetic state within 
the canted antiferromagnetic medium. It appears canted, 
the spins inside the droplets being deviated from the per- 
fect AF structure by an angle of «18 degrees (9 av «10 
degrees). The temperature variation of the intensity be- 
low Tjv, must be related to the evolution of the contrast 
Am between the two magnetic phases, which varies with 
the mean magnetization. The almost temperature inde- 
pendent evolution of q m suggests that the charge carrier 
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FIG. 3. (a) Energy scans fitted by energy-lorentzians for 
q|| [00£] at 1 6K. From top to bottom, £ varies from 2.05 to 2.5. 
(b) Dispersion curves for x=0.05 (filled circle) and x=0.08 
(filled square). The continuous line at upper energy corre- 
sponds to the spin wave branch observed in pure LaMnC>3. 



segregation persists at Tjv- Above T^r, where the spa- 
tial distribution described above cannot be observed, the 
information arises from the spin dynamics only, which in- 
dicates a ferromagnetic correlation length typical of these 
inhomogeneities |T^ ]. 

We compare now the characteristics of the static mag- 
netic inhomogeneities with those of the low energy spin 
wave branch which appears by doping the system, con- 
sidering both its q-dependent intensity and its dispersion 
curve. In Fig. 3 (a) we have reported some energy spec- 
tra obtained as a function of q//[001]. Similar results 
are observed for q//[110], indicating an isotropic and 
weakly dispersed spin wave branch. As previously shown 
for x=0.05 fig ], the intensity of this spin wave branch 
strongly decreases with q. Fitting the dynamical form 
factor (or energy integrated intensity) with a lorentzian 
function, leads to determine a ferromagnetic correlation 
length f with £ « 7.5 - 8l for both [110] and [001]. This 
reveals a tight connection between this low energy spin 
wave branch, characteristic of a low and isotropic cou- 
pling, with the static droplets described above. Among 
other characteristic features, these spin excitations keep 
a propagative character and are still well-defined even at 
small q, i. e. over distances larger than the droplet size. 
These excitations may therefore correspond to eigen- 
states of the whole spin system. Within the above in- 
homogeneous picture derived from the static study, this 
means that the droplet magnetic states are coupled to- 
gether through the matrix. In Fig 3-b, the dispersion 
curve is shown along [001] together with that of the other 
spin wave branch. The two spin wave branches observed 
for x=0.05 are also reported for comparison. When con- 
sidering the dispersion curve of the low energy branch, a 
fit using w=Dq 2 +cjo as an approximation at small q, indi- 
cates an increase of the stiffness constant D with x, and a 
decrease of the gap at t=(002). Such an evolution of the 



constant D with the number of carriers can be expected 
at larger doping, within the true metallic and ferromag- 
netic state pi] . However, within the present inhomo- 
geneous picture, we observe a very peculiar evolution of 
the dispersion curves with x. The corresponding two spin 
wave branches keep separated by an energy gap, constant 
with x (hatched area in Fig 3-b). Such a pecularity also 
holds at x=0.1. Therefore, the separation in the direct 
space into two magnetic states, induced by an electronic 
phase segregation is associated with a separation in the 
energy space into two spin dynamics, characteristic of 
two types of magnetic coupling. 

In conclusion, the very unusual features found for the 
static and dynamical spin correlations provide new in- 
sights in the physics of the electronic phase separation. 
A liquid-like spatial distribution of magnetic polarons has 
been observed for the first time. The picture of hole-rich 
droplets agrees with predictions of the s-d model in a 
well-known limit ||||J|- Whether such inhomogeneities 
exist at x<x c where the compound becomes metallic and 
ferromagnetic, requires more studies. 

One author (M. H.) is very indebted to B. Hennion 
and D. Khomskii for fruitful discussions and G. Coddens 
for his critical reading of the manuscript. 
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